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ABSTRACT 

We present new Chandra and XMM-Newton observations of a sample of eight radio-quiet 
7-ray pulsars detected by the Fermi Large Area Telescope. For all eight pulsars we identify the 
X-ray counterpart, based on the X-ray source localization and the best position obtained from 7- 
ray pulsar timing. For PSR J2030-I-4415 we found evidence for a ^ 10"-long pulsar wind nebula. 


Our new results consolidate the work from Marelli et al. (20111 and confirm that, on average, the 


7-ray-to-X-ray flux ratios (F.y/Fx) of radio-quiet pulsars are higher than for the radio-loud ones. 
Furthermore, while the F.y/Fx distribution features a single peak for the radio-quiet pulsars, the 
distribution is more dispersed for the radio-loud ones, possibly showing two peaks. We discuss 
possible implications of these different distributions based on current models for pulsar X-ray 


emission. 


Subject headings: Stars: neutron — Pulsars: general — X-rays: stars — gamma rays: stars 


1. Introduction 

The launch of the Fermi 7-ray Space Telescope 
in 2008 marked a revolution in pulsar 7-ray as- 


ond and young-to-middle-aged rotation-powered 
7-ray pulsars, with 117 listed in the Second Fermi 
LAT Catalog of 7-ray pulsars (2PC, |Abdo et~^ 
2013[ ). This number has already risen to more 
than 16(Q About 30% of the LAT pulsars are 
radio-quiet. In the 2PC, radio-quiet (RQ) pulsars 
are defined as those that have not been detected in 
radio down to a flux density limit of S1400 = 30/rJy 


^ http: //tinyurl.com / fermipulsars 


tronomy (for a recent review see Caraveo 20141. 


The Large Area Telescope (LAT, Atwood et al. 


20091 onboard Fermi has detected tens of millisec- 


at 1400 MHz, whereas radio-loud (RL) pulsars 
are those detected above this limit. Some pulsars 
that have been detected in radio with a flux den¬ 
sity below Si4oo = 30^Jy are defined radio faint 
(RF). Assuming that radio emission is produced 
close to the magnetic poles, this large fraction of 
RQ 7-ray pulsars suggests that the 7-ray beam is 
broader and at a larger angle from the magnetic 
poles with respect to the radio beam (seelCaraveo 


2014, and references therein), thus making its de¬ 
tection less sensitive to geometric effects. This 
piece of evidence has been crucial to verify the 
predictions of pulsar magnetospheric models, such 


as the outer gap and the polar gap ones (Cheng 


et al. 1986 Harding&Muslimov 20041. Multi¬ 


wavelength studies of 7-ray pulsars are key to 
map the geometry of the different emission regions 
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in the pulsar magnetosphere, investigate possible 
connections between different emission processes, 
and study their efficiencies as a function of en¬ 
ergy. For instance, [Marelli et al. (2011); Marelli 


( 2012 ) showed that the distribution of the 7 -to- 
X-ray flux ratio Fly/Fx is different for RL and RQ 
pulsars, being narrower and peaking at higher val¬ 
ues for the latter. What is behind the difference 
in the relative 7 and X-ray efficiencies of these 
two classes of 7 -ray pulsars is not clear yet. While 
some clues are evident from the study of a handful 


of RQ pulsars detected in X-rays (see e.g. Marelli 
et al.|[2M3 2014a|b ), studying a larger sample al¬ 


lows one to explore differences across the parame¬ 
ter space. Furthermore, the detection of RQ 7 -ray 
pulsars in X-rays (and in the optical) can provide 
information, for instance, on the pulsar dynamics 
and distance, traditionally obtained in the radio 
band. In fact, for RL pulsars the measurement of 
the dispersion of the pulses at different radio fre¬ 
quencies allows one to estimate the free electron 
column density, from which the distance to the 
pulsar is obtained from a gas distribution model 
(Cordes & Lazio 2002). Alternatively, direct dis¬ 
tance measurements for RL pulsars are obtained 
from the radio pulsar parallaxes. The measure¬ 
ment of the proper motion and distance of the RQ 
LAT pulsar PSR J0357-I-0352 with Chandra and 
XMM-Newton ( De Luca et al.|p013 Marelli et al. 


2013) is a spectacular example of distance estima¬ 


tion is possible also for RQ pulsars and using a 
different method. 


The X-ray satellites that have yielded the high¬ 
est pulsar detection rates have traditionally been 
Chandra and XMM-Newton due to their spatial, 
timing and spectral resolutions, along with their 
low background rates. Here, we present the re¬ 
sults of new follow-up observations of a sample 
of RQ LAT pulsars with no previous X-ray detec¬ 
tions. Our observations are summarized in Section 
2 , while our data analysis and results are described 
in Section 3. In Section 4 we compare the X and 
7 -ray properties of all RQ LAT pulsars, making 
some considerations based on current models of 
7 -ray and X-ray emission from pulsars. 


2. Target selection and observation de¬ 
scription 


We selected our target pulsars for Chandra and 
XMM-Newton observations among those RQ LAT 
pulsars that have no, or uncertain, detection in 
X-rays (see Table 5 of 2PC). In order to predict 
the non-thermal X-ray flux of these pulsars, we 
relied on the observed Fy/Fx distribution of the 
RQ pulsars family, favoring pulsars with relatively 
high 7 -ray flux and small pseudo-distance Dy, as 
inferred from the comparison between flux and 
luminosity of LAT pulsars with known distance 
( Saz Parkinson et al.||2010 ). Our sample includes 
eight LAT pulsars, spanning three orders of mag¬ 
nitude in characteristic ages (tc = F/2F, where P 
is the neutron star spin period and P is derivative) 
and two orders of magnitude in spindown energies 
(F = Aid^lP/P^ , where the moment of inertia I is 
assumed to be 10^® g cm“^). Their 7 -ray timing 
coordinates, which we used as a reference for the 
X-ray counterpart identification, and their timing 
parameters are summarized in Table 1, together 
with the values of the pseudo-distance Dy, as de¬ 


fined in Saz Parkinson et al. (2010). We took the 


pulsar timing parameters and the pseudo-distance 
values from the compilations in Abdo et al. (2013), 
whereas we took the 7 -ray timing positions and er¬ 
rors from the most recent compilation of the LAT 
7 -ray Pulsar Timing Models This compilation 
reports the results of the analysis of five years of 
Fermi data by using a new, advanced timing anal¬ 
ysis pipeline ( Kerr et al.jin preparation ) that han¬ 
dles glitch detection and fitting in the presence of 
timing noise much better than previous codes. We 
note that before this work, the positional errors 
based on 7 -ray timing were correctly calculated 
(e.g. taking into account timing noise and sys¬ 
tematic errors) only for a small number of bright 
pulsars ( Ray et ^|2011 see e.g.). 

Swift already observed all the 7 -ray pulsars in 
our sample in short (« 5-10 ks) snapshot obser¬ 
vations. Through the re-analyses of these data, 
we found no X-ray source that could be position- 
ally associated with the pulsar (for a more detailed 
discussion see Marelli ||2012 ). Three of these pul¬ 
sars were also re-observed by Suzakw. two of them 
(PSR J1429-5911 and J1957-b5033) remained un- 


^https://confluence.slac.stanford.edu/display/ 
GLAMCOG/LAT+Gamma-ray+Pulsar+Timing+Models 
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detected, also owing to the short exposure times, 
whereas a possible, marginal X-ray detection was 
obtained for PSR J1838-0537. 

Here, we report on five new 25 ks Chandra 
observations of PSRs J1429-5911, J1957-h5033, 
J2028-h3332, 32030-^4415, and 32139-^4716 (ob¬ 
servation ids 14825, 14828, 14826, 14827 and 
14829, respectively) and 15 ks new Chandra ob¬ 
servations of PSRs 30734—1559 and 31846-1-0919 
(observation ids 13792 and 13793, respectively). 
We observed all the pulsars with the Advanced 
CCD Imaging Spectrometer (ACIS). In order to 
detect the highest number of counts from our tar¬ 
gets we obtained the ACIS-S observations in the 
VFAINT mode, with the targets positions placed 
on the back-illuminated ACIS S3 chip. We ob¬ 
served PSR 31838—0537 with XMM-Newton for 
44 ks (obs. id 0720750201) using The EPIC (Eu¬ 
ropean Photon Imaging Camera) PN and the two 
Metal Oxide Semi-conductor (MOS) cameras in 
the Pull Erame mode with medium optical filters, 
due to the presence of moderately bright stars 
within their field-of-view (POV). All the obser¬ 
vations were carried out between 2012 September 
16 and 2014 April 15. 

3. Data analysis and results 

Eor the Chandra data analysis we used the 
Chandra Interactive Analysis of Observation 
(CIAO) software version 4.3. We re-calibrated our 
event data by using the chandra_repro script and 
selected events in the 0.3-10 keV energy range. 
Depending on the total number of counts, we ran 
the wavdetect detection tool in different energy 
ranges, taking into account the exposure maps. 
Eor each pulsar, we found a single X-ray source 
within (or close to) the 95% 7 -ray timing error 
box, with a detection significance above 5a and 
source counts varying from 12 (PSR 31846-1-0919) 
to 123 (PSR 31957-1-5033). Eor each X-ray source 
we performed a spectral analysis. We extracted 
the source counts from a 2 " circle radius around 
the X-ray position and the background counts 
from a surrounding annulus, whose radius we 
chose on a case by case basis in order to avoid 
contamination from serendipitous nearby sources. 
We used the CIAO tool specextract to create 
the spectra, response matrix and effective area 
files and analyzed the spectra with XSPEC (ver¬ 


sion 12.8.1). We used the C-statistic ( Cash[19TO l, 
based on the application of the likelihood ratio and 
recommended for cases with low statistics and low 


background. As done in Marelli et al. (20131, we 


performed a brightness profile analysis on all X- 
ray sources, revealing no evidence for extended 
emission, with the exception of PSR 32030-1-4415 
(see Section 3.7). 

Eor the analysis of the 44-ks XMM observation 
of PSR 31838—0537, we used the XMM-Newton 
Science Analysis Software (SAS) vl3.0. We per¬ 
formed a standard analysis of high particle back¬ 
ground (following De Luca et al.|[2005 1 and cross¬ 
checked the results with the SAS tool bkgoptrate 
(also used for the 3XMM source catalog, obtain¬ 
ing very good agreement. After the subtraction 
of bad time intervals affected by soft proton X- 
ray flares, we obtained a net exposure time of 23 
ks. We selected 0-4 pattern events for PN and 
0-12 for the MOS detectors in the 0.3-10 keV en¬ 
ergy range. Prom the cleaned events, we ran the 
source detection according to two different meth¬ 
ods: we used the SAS edetect-chain and the 
CIAO wavdetect tools. Both of them detected 
a single X-ray source within few arcsecs from the 
pulsar 7 -ray timing position. The source appears 
point-like, with its brightness profile consistent 
with the detector’s Point Spread Punction. In or¬ 
der to maximize the source signal-to-noise (S/N), 
we extracted spectral counts from a 25" radius cir¬ 
cular region around the computed X-ray position 
and the background from a source-free, nearby, 
40" radius circle on the same CCD. We generated 
ad hoc response matrices and effective area files 
using the SAS tools rmfgen and arfgen. To in¬ 
crease the statistics, we added the two MOS spec¬ 
tra by using the HeaSoft tool mathpha and the two 
response matrices and effective area with addrmf 
and addarf. After background subtraction, we 
extracted 200 and 176 net source counts from the 
PN and the two MOS, respectively. Owing to the 
high background (64% and 43% of the extracted 
counts from the PN and MOS, respectively), we 
used the statistic for the spectral fit. 

Due to the low statistics of our X-ray sources, 
we assumed the pulsar pseudo-distance to nor¬ 
malize the value of the integrated Galactic Nh in 


^http://xmmssc-www.star.le.ac.uk/Catalogue/xcat_ 
public_3XMM-DR4.html 
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the pulsar direction, computed according to the re¬ 


calibration (Schlafly&Finkbeiner 20111 of the ex¬ 


tinction maps from Schlegel et al. (19981. Four of 


our eight sources have fewer than 20 net counts, 
which are not adequate to perform detailed spec¬ 
tral fits. Therefore, for these sources we fixed the 
power-law (PL) photon index and blackbody (BB) 
temperature to representative values of 2 and 200 
eV, respectively. Following, e.g., [Marelli et al.| 
(20111, these are the fitted average of the values 


measured for pulsars detected in the X rays. Er¬ 
rors in the spectral parameters are reported at the 
90% confidence level (c.L). 

The errors on the X-ray positions reported in 
the following sub-sections are purely statistical, at 
a 3a confidence. To them, we have to add the 
90% c.l. systematical errors associated with the 
absolute accuracy of the satellite aspect solution, 
which are Ol'S and 1"5 per coordinate for Chan- 
drc0 and XMAi-Newtoi^ respectively. According 
to the logN-logS distribution of Chandra Galac¬ 
tic sources ( Novara et al.||20d9 ), we can estimate 
the probability of a chance detection of an X-ray 
source within a representative LAT y-ray timing 
error box (1" see Figure [^, with X-ray flux 
similar or greater than the measured ones, to be 
about 3x10“^. Thus, based on positional coin¬ 
cidence, we consider our identifications to be se¬ 
cure. Figure shows the positions and errors of 
y-ray pulsars and the associated X-ray counter¬ 
parts. The computed X-ray spectral parameters, 
unabsorbed X-ray fluxes, and y-to-X-ray flux ra¬ 
tios of our eight pulsars are summarized in Table 
2 . 

3.1. PSR J0734-1559 

PSR J0734—1559 (P=156 ms) was identified 
as a middle-aged (0.2 Myr) y-ray pulsar during 
a blind search for pulsations from the unidenti- 


son 


tied LAT source IFGL J0734.7-1557 (Saz Parkin- 

12 ^. 


We detected its X-ray counterpart at 
ax = 07*'34“45!7 (±(y.'15); Sx = -15°59'19"8 
(±Cf.^26), consistent with the pulsar y-ray timing 
coordinates. The X-ray source is detected with a 
significance of 8.3 a, as computed by wavdetect. 


estimate a Nh = 2x 10^^ cm“^. Keeping these val¬ 
ues fixed, we fitted the X-ray spectrum of the pul¬ 
sar with either a single power law (PL) or black- 
body (BB) model. Due to the low statistics, we 
also fixed the PL photon index to 2 and the BB 
temperature to 200 eV. The best-fit with a BB 
gives a radius of the emitting region on the neu¬ 
tron star surface of m, computed for the 

assumed pseudo-distance D.^. For the PL model, 
the unabsorbed X-ray flux in the 0.3-10 keV en¬ 
ergy range is = 1.5±0.6x 10“^^ erg cm“^ s“^, 
which gives a y-to-X-ray flux ratio /Fx ^ 3700. 
The BB model, on the other hand, gives an X-ray 
flux Fx® = 1.2 ±0.4 X 10“^^ erg cm“^ s“^ and an 
F-y/Fx > 4700, assuming F®® as an upper limit 
on the non-thermal X-ray flux from the pulsar. 

3.2. PSR J1429-5911 

The y-ray pulsar PSR J1429—5911 (P=115 ms) 
was one of the very first discovered by applying 


the blind-search technique (Saz Parkinson et al. 


2010). We detected its X-ray counterpart at ax = 
TW9'"58!5 (±1"03); 5x = -59°lT36"2 (±0"45), 
with a significance of 8.5a (24 net source counts). 
As we did in the previous section, we assumed 
the pseudo-distance D.y=1.7 kpc to estimate an 
Nh = 3 X 10^^ cm“^ and we kept it fixed in the X- 
ray spectral fit. In this way, a fit with a PL gives a 
photon index Fx = —0.1±0.7 and an unabsorbed 


X-ray flux F®® = (3.3 ± 1.8) x IQ-^^ erg cm-^ 
s“^. This corresponds to Fy/Fx ^ 2400. A fit to 
the data with a single BB model did not yield an 
acceptable temperature (T>0.8 keV). 

3.3. PSR J1838-0537 

The y-ray pulsar (P=145 ms) PSR J1838—0537 
was discovered through a blind search for pul¬ 
sations of the unassociated LAT source 2FGL 
J1839.0—0539 (iPletsch et al.||2012a|). The pul¬ 


sar is the youngest (5 kyr) and most energetic 
(F ~ 5.9 X 10^® erg cm“^ s“^) in our sample. 
PSR J1838—0537 is spatially coincident with the 
TeV source HESS J1841—055, hence possibly as¬ 
sociated with a PWN detected at very high ener¬ 
gies. A 41.1 ks Suzaku observation revealed a can- 


with 19 net counts. As described above, we as- 

didate X-ray counterpart ( 

Pletsch et al. 

2012al 

sumed the pulsar pseudo-distance D.y=1.3 kpc to 

However, the non-negligible probability o 

' a find 


http://cxc.harvard.edu/cal/ASPECT/celmon/ 
^calibration technical note XMM-SOC-CAL-TN-0018 


ing a spurious source within the large Suzaku er¬ 
ror circle (« 19" radius Uchiyama et al.|[2008l, as 


well as the low significance of detection (~ 3a) 
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made the identification of the X-ray source with 
the pulsar uncertain. Owing to the better angu¬ 
lar resolution of XMM-Newton we detected the 
pulsar X-ray counterpart at ax = 18*'38™56!2 
(±2"65); Sx = -05°37'04"5 (±2'.'69), with a sig- 
nihcance of 20.7 a (~ 270 counts; pn-|-MOS). The 
X-ray spectrum is best fitted by an absorbed PL 
with a photon index Lx = 0 . 8 i'i|!) g and a fitted 
Nh = X 10^^ cm-2 (90% c.L). The un¬ 

absorbed X-ray flux is = (7.2 ± 0.9) x 10“^^ 
erg cm“^ s“^, which gives a y-to-X-ray flux ratio 
^ 7 /Fx 2600. A fit with a BB spectrum did not 

yield an acceptable temperature (T>1.5 keV). 


3.4. PSR J1846+0919 


PSR J1846-f0919 (P= 225 ms) was identi- 
fied as a 7 -ray pulsar by |Saz Parkinson et al. 
(20101. We detected the X-ray counterpart at 
ax = 18*'46“25!8 (±0'.'41); ,5x = 09°19'49"8 
(TOi'dS) with a signihcance of 5.4cr (12 net source 
counts). As usual, from the pseudo-distance we 
estimated an Nh = 2 x 10 ^^ cm“^ and kept it 
fixed in the X-ray spectral fit. A fit with a PL 
with fixed photon index (Lx = 2 ) gives an un¬ 
absorbed X-ray flux = (8.7 ± 4.8) x lO-^^ 
erg cm“^ s“^, that corresponds to F^/Fx 2800. 
Owing to the low statistics, the X-ray spectrum 
of this source can be fitted also by a BB. From 
the pseudo-distance of 1.4 kpc we obtained an 
emitting radius of 75^i[g m and the unabsorbed 
flux is Ff ^ = (4.5 ± 3.3) X lO-^^ erg cm-^ s"!. 
Assuming, as before, as an upper limit to the 
non-thermal X-ray flux of the pulsar, the lower 
limit on F^/Fx would be 5300. 

3.5. PSR J1957+5033 

Like PSR J1429-5911 and J1846-f0919, PSR 
J1957-I-5033 (P=374 ms) is one of the very first 
pulsars discovered by a blind search by |Saz Parkin-] 
(|2010D. We detected the X-ray coun- 
= 19*'57“38M (±(y.'20); 5x = 
The pulsar X-ray counter- 


son et al. 


terpart at ax 
50°33'20'.'8 (±0'.'17). 
part is detected with a signihcance of 29 .Oct (123 
net source counts). Owing to the adequate statis¬ 
tics, we left the absorption column as a free pa¬ 
rameter in our spectral hts, as in the case of PSR 
J1838—0537 (Section 3.3). The X-ray spectrum is 
best ht by a PL, which gives an Nh < 2.5 x 10^° 
cm“^, a photon index Px = 2.1 ± 0.3, and an un¬ 


absorbed flux Fjf^ = (3.0±0.5) X 10“^"* erg cm“^ 
s“^. This corresponds to a quite low F^/Fx ~ 
870, indeed the lowest among the pulsars in our 
sample. A ht with a single BB component is not 
acceptable (null hypotesis probability of 4x 10“®). 


3.6. PSR J2028±3332 


This pulsar (P=176 ms) is one of the nine 7 -ray 
pulsars discovered by Pletsch et al. (2012b I us¬ 


ing a novel blind search technique. We identihed 
the PSR J2028±3332 X-ray counterpart at ax = 
20*'28“19’!8 (±0"45); (5x = 33°32'04"1 (±0"23) 
with a signihcance of 5.1 ct (15 net source counts). 
As we did in previous cases, we used the pseudo¬ 
distance Dg,=0.9 kpc to derive an Nh = 2 x 10^^ 
cm“^. The ht with a PL with hxed photon in¬ 
dex (Lx = 2 ) gives an unabsorbed X-ray hux 

and 


= (5.3 ± 3.3) X 10 erg cm 


pPL 

Fj/Px ^ 10900, the highest among the pulsars in 
our sample. We also tried a ht with a BB, which 
gives an emitting radius of only 34lg^ m, an X-ray 
hux F®® = (2.2 ± 1.5) X 10“^® erg cm“^ s“^, and 
an Fg,/Fx > 26400. 


3.7. PSR J2030±4415 

The 7 -ray pulsar J2030±4415 (P=227 ms) was 
also discovered 


by 


(2012b). The 


Pletsch et al. 
pulsar X-ray counterpart, at ax = 20*'30™5T4 
(±0'.'15); 6x = 44°15'38"8 (±0'.'16) is detected 
with a signihcance of 15 ct (54 net source counts). 
The 0.3-10 keV X-ray spectrum is described by a 
PL with with photon index Px.psr = for 

an Nh = x 10 


20 , 


cm 

sorbed X-ray hux Fx.psr = ( 2.1 ± 0 . 6 ) x 10 


, which gives an unab- 

— 14 

erg 

cm“^ s“^. The low value of the absorption column 
(Nh ^ 3 X 10^^ cm“^) with respect to the Galactic 
Nh in the pulsar direction (^ 10 ^^ cm“^) agrees 
with the pseudo-distance Dg,=800 pc. The 7 -to- 
X-ray hux ratio for PSR J2030±4415 is F^/Fx ~ 
2800. A ht with a single BB component is not ac¬ 
ceptable (null hypothesis probability of 6 x 10 “®). 
We found a ^ 10"-long extended emission around 
the pulsar, elongated North to South. To search 
for extended emission, we applied the CIAO tool 
vtpdetect on an event list purged from the point¬ 
like sources we found in the held using wavdetect. 
This tool, based on the Voronoi Tessellation and 


Percolation (VTP, see e.g. Boschin 2002), is partic¬ 
ularly indicated for the search of extended sources. 
This resulted in a false source probability (i.e. the 
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probability that the detection is associated to a 
real source) for PSR J2030+4415 nebula of 10“^® 
and a best-fitting elliptical region with a major 
axis of ^ 10". Assuming the pseudo-distance of 
800 pc, the ’’tail” would have a physical dimen¬ 
sion of - 0.07 pc. Since PSR J2030-h4415 is RQ 
and Chandra has observed it only once, we have 
no proper motion information yet, thus we can¬ 
not say whether the X-ray nebula is aligned with 
the motion of the pulsar, although this would be 
the most likely interpretation. We extracted the 
net counts from the PWN from an elliptical re¬ 
gion of 10" semi-major axis after masking the pul¬ 
sar. We computed the background from nearby, 
source-free elliptical regions. The tail spectrum is 
well described by a PL, with Px.pwn = 
with the unabsorbed X-ray flux of Px,pwn = 
(4.0 ± 1.6) X 10“^"^ erg cm“^ s“^, about twice as 
large as the pulsar. 


3.8. PSR J2139+4716 


This 7 -ray pulsar (P=282 ms; Pletsch et al. 
2012b) is the oldest (2.5 Myr) and least energetic 
{E ~ 0.3 X 10^“^ erg cm“^ s“^) in our sample. We 
detected the pulsar X-ray counterpart at ax = 
21*'39“56!0 (±01'44); 5x = 47°16'13"0 (±0"52) 
with a significance of 5.5 (t (16 net source counts). 
As usual, from the pseudo-distance D.y=0.8 kpc we 
derived Nh = 10^^ cm“^. By keeping it fixed, a fit 
with a PL with fixed photon index (Px = 2) gives 
an unabsorbed X-ray flux = 4.7± 2.5 x 10“^® 

erg cm“^ s“^, and F^/Fx 4900. The fit with 
a BB gives a very small emitting radius of only 
33 ± 9 m, a X-ray flux = (2.7 ± 1.4) x lO-^^ 
erg cm“^ s“^, and a F^/Fx> 8500. 

3.9. Optical and infrared observations 

None of the pulsars in our sample have been 
observed in the optical or infrared (IR, 2PC). 
For completeness, we also scanned optical/IR 
data from public imaging surveys. For PSR 
J1429—5911 we found serendipitous IR observa¬ 
tions from the VVV (VISTA Variables in the Via 


Lactea; Emerson et al. 2006) survey, carried out 


at the ESO’s Cerro Paranal Observatory (Chile) 
with the 4.1m Visible and Infrared Survey Tele¬ 
scope for Astronomy (VISTA) and an IR camera 
(VIRCAM; [Dalton et aljM^ . The fields of PSR 
J1838—0537 and PSR J1846±0919 were serendipi- 
tously observed in the UKIDSS (UK Infrared Deep 


Sky Survey; Lawrence et al.||2007 ), carried out at 
the Mauna Kea Observatory (Hawaii) with the 3.8 
m UK Infrared Telescope (UKIRT) and the Wide 
Eield Camera (WFCAM; Casali et al.|2007 ). In all 
these cases, no objects are detected at the Chan¬ 
dra or XMM-Newton positions. The derived up¬ 
per limits are orders or magnitudes above the flux 
levels expected for the pulsar’s age and spin-down 
power by scaling for the luminosity and distance 


of other pulsars detected in the IR (e.g. Mignani 


et al. 2012). Thus, these results are mainly to be 


considered as a reference for future deeper follow¬ 
up observations. Our 3 (t limits are Z>21; Y>21; 

20; H>19; K > 19 (PSR J1429-5911), J~17, 
H-19, K-18 (PSR J1838-0537), and J- 20.5, 
H- 19.4, K~ 18.6 (PSR J1846±0919). 


4. Discussion and conclusions 

Our Chandra and XMM-Newton observations 
yielded the identification of the X-ray counterparts 
of eight RQ LAT pulsars. Thus, the number of RQ 
LAT pulsars detected in the X-rays amounts now 
to 28, with only 11 still missing an X-ray coun¬ 
terpart. This number has to be compared with 
the 30 RL LAT pulsars (out of 42) detected in the 
X-rays. Our enlarged X-ray database is important 
for a better understanding of the multi-wavelength 
emission properties of LAT pulsars. In particu¬ 
lar, focussing on the differences between young- 
to-middle-aged RQ and RL pulsars, |Marelli et al.| 
(2011) found that, for a given y-ray flux, the for¬ 


mer are intrinsically fainter in the X-ray band than 
the latter. Thus, the distance-independent F.y/Fx 
turned out to be a useful parameter to investigate 
the differences between the X- and y-ray emission 
properties of the two pulsar families. Eor consis¬ 
tency, we followed the same approach as in Marelli 
et al. (2011) and we used our much enlarged sam¬ 
ple of X-ray-detected RQ pulsars to improve their 
analysis. 

We note that millisecond pulsars (MS) were not 
included in the analysis of Marelli et al. ( 2011| ), 
owing to the fundamental differences between 
their magnetospheres and those of young/middle- 
aged pulsars (e.g. Johnson et al. 2014). In fact, ow¬ 


ing to their shorter periods, MS pulsars are though 
to have a smaller magnetosphere than young pul¬ 
sars, where the radius of the last closed field line is 
defined by the so-called light-cylinder radius Rlc 
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= P cj^'K. Moreover, given their much longer and 
more complicated evolutionary history, the mag¬ 
netic field topology in MS pulsars may be more 
complex than those of younger pulsars. This may 
in part explain why their pulse profiles are more 
intricate and diverse. Lastly, MS pulsars gener¬ 
ally seem to be more massive than their younger 
counterparts. See Venter et al. (2014) for a more 
detailed discussion. 

In addition to those obtained from our obser¬ 
vations, we collected 7 - and X-ray best-fit spec¬ 
tra from the 2PC and papers published after¬ 
wards (PSR J1357-6429, |Chang et al.| (I20T2|; 
J1741-2054, IMarelli et al.l (120144; J1813-1246, 


Marelli et al. (|2014bD; J0357-f3205, 

Marelli et 

al. (2013); J2055-b2539, 

Marelli et al. 

(in prepa- 


Log(F.y/Fx) values for all pulsars, both RL and 
RQ, detected in the X-rays and with a non- 
thermal X-ray spectrum. Improving on the 
method used in Marelli et al. (2011); Marelli 


( 2012 ), each pulsar has been represented with an 
asymmetric parabola, in order to account for the 
asymmetric errors on the X- and 7 -ray fluxes. We 
separately highlighted in Figurej^the Log(F.y/Fjf) 
values of the eight pulsars with the newly discov¬ 
ered X-ray counterparts (Table 2), assuming a 
non-thermal X-ray spectrum for all of them. We 
fitted both RQ and RL pulsar distributions with 
single and double gaussian models, taking into 
account Poissonian errors. 


We confirmed that RQ pulsars are intrinsically 
fainter in the X rays than the RL ones, with an 
average Log(F.y/Fx) of 3.38±0.10 (Icr error): at 
mean, radio-loud pulsars have F.y/Fx values an 
order of magnitude lower (see later). The dis¬ 
tribution of the F-y/Fx values for RQ pulsars is 
well-fitted by a gaussian, with a null hypothe¬ 
sis probability of 0.62, three degrees of freedom 
(dof), featuring a sharp peak (standard deviation 
of 0.43±0.09). This indicates very similar X- and 
7 -ray emissions among the members of this family. 
The F.y,/Fx values computed for the eight RQ pul¬ 
sars with newly discovered counterparts (Table 2) 
are distributed around the peak, significantly in¬ 
creasing our statistics. We note that for four of 
the eight pulsars the X-ray spectrum could also 
be fitted by a single BB model, which would re¬ 
sult in an even higher value of the corresponding 
F^/Fx. 


As noted from Figure]^ the F.y/Fx values of 
RL pulsars are characterized by a more struc¬ 
tured distribution. While a single gaussian could 
fit the distribution, with an average Log(F.y,/Fx) 
of 2.24±0.32, the null hypothesis probability re¬ 
sults quite low (0.05, eight dof) and the distri¬ 
bution is much wider, with a standard deviation 
of 0.72±0.11. By fitting the distribution with two 
gaussians, we obtain two peaks at I.8I±0.1I (stan¬ 
dard deviation of 0.29±0.09) and 3.20±0.14 (stan¬ 
dard deviation of 0.34±0.05), with the better null 
hypothesis probability of 0.38 (five dof). An f-test 
( Bevington||I969 ) shows that the probability for a 
chance improvement is 0.03, not enough to exclude 
the single-gaussian fit. 


Interestingly enough, the second peak of the 
F-y/Fx distribution of RL pulsars would overlap 
the peak of the corresponding distribution for the 
RQ pulsars. If statistically confirmed, this might 
suggest that the separation between RL and RQ 
pulsars in the F.y,/Fx space might not be as clear 
as it was originally thought. At least partially, 
the overlap between the two peaks might be at¬ 
tributed to the subtle, and somehow arbitrary, dis¬ 
tinction between radio-faint and radio-quiet pul¬ 
sars. Indeed, some RL pulsars with the highest 
F^/Fx values (e.g., J1741-2054 and J1907-f0602, 
see Figure 3 in the 2PC) are radio-faint, whereas 
some RQ pulsars with the lowest F.y/Fx values, 
e.g. J18I3—1246, are very distant and absorbed. 
Thus, they would be undetected if they had a ra¬ 
dio luminosity comparable to those of J1741—2054 
and J1907-I-0602. In any case, we found no clear 
correlation between the radio flux and the X- and 
7 -ray fluxes, which might have suggested that, e.g. 
radio-faint pulsars have higher F.y/Fx values than 
the radio-bright ones. 


In a similar way, we built the histogram of the 
photon indices ratio, r.y/rx, for RL and RQ pul¬ 
sars families (Figure]^. In both cases, the dis¬ 
tributions are well fitted by a gaussian (null hy¬ 
pothesis probabilities of 0.16 and 0.62, 2 dof for 
RL and RQ pulsars, respectively). Interestingly 
enough, the peaks of both distributions occur at 
very similar values of the F..y/Fx ratio, -0.09±0.03 
and -0.14±0.08 for the RL and RQ pulsars, respec¬ 
tively. This suggests that the peculiar distribution 
seen in the F.y/Fx histogram (Figure]^ is not pro¬ 
duced by an intrinsic difference between the spec¬ 
tral slopes of the two families. In the following. 
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we explore some implications of the observed bi- 
modal F^/Fjf distribution for RQ and RL pulsars 
on pulsar emission models. 

Pulsar magnetospheric radiation is highly 
anisotropic and a complex antenna pattern (i.e., 
the direction dependence of the emitting power) 
results as the neutron star rotates and its beam 
sweeps the sky. In a specific energy range, dif¬ 
ferent emission models expect different antenna 
patterns - the emission mechanism, geometry and 
luminosity depend on the magnetic field configu¬ 
ration, inclination angle and intensity, as well as 
on the pulsar period, but the physics ruling such 
mechanisms is not yet understood. For a complete 


discussion see e.g. Bai&Spitkovski (2010). Con¬ 


straining the multi-wavelength antenna patterns 
as a function of different pulsar properties would 
yield crucial clues to understand pulsar magneto¬ 
spheres. 

The observed flux in a given energy range is a 
phase-averaged cut through the antenna pattern 
for the Earth line of sight. A beaming factor /n 
is usually defined as the ratio of the observed flux 


to the average flux over 47 r sr (e.g. Watters et al. 


2009 1 . In the radio range, the antenna pattern is 


usually described with the ” cone plus core” heuris¬ 


tic model (e.g. Harding 2007). According to this 
model, then, the radio emission is centered on the 
magnetic axis (in a co-rotating frame) and radio 
loudness (or, quietness) is merely a fortuitous ef¬ 
fect, set by our line of sight intercepting (or, miss¬ 
ing) the radio antenna pattern of a pulsar, as a 
result of a small (or, high) magnetic impact an¬ 
gle. In other words, RQ pulsars are simply those 
that are seen under a viewing angle substantially 
different from their magnetic inclination angle - 
the beaming factor for such lines of sight being 
essentially zero. 

In the 7 -ray range, recent works based on data 


collected by the Fermi-LMJ: (see e.g. Bai&Spitkovski 


2010 Pierbattista et al. 2014) showed very com¬ 


plex antenna patterns, best described by outer 
magnetospheric models such as the Outer Gap 


(Romani 1996) or the Two-pole Caustic model 
(Muslimov & Harding 2004| ). According to such 
models, RQ pulsars should not have, on average, 
a different 7 -ray beaming factor with respect to 


RL ones (e.g. Pierbattista et al. (2014) - indeed. 


beaming factors for RQ pulsars are expected to be 
more dispersed than for RL pulsars). On the other 


hand, although distance estimates for RQ pulsars 
are highly uncertain, there is no evidence for a 
larger luminosity (nor for a larger spin-d own con- 
version efficiency n.y =L.y/A, as defined in 


Abdo et 


al. (2010|) with respect to the RL ones, although 


some selection bias, related to highly pulsed 7 -ray 
signal and little timing noise easing blind period¬ 
icity searches in 7 -ray data, is certainly affecting 
our view of RQ pulsars. 

In the above picture, it is tempting to link the 
different ¥^/Yx distributions for RQ and RL pul¬ 
sars to differences in their X-ray emission proper¬ 
ties, which could allow us to set first constraints on 
the poorly known X-ray antenna pattern. For the 
RQ pulsars (seen with large magnetic impact an¬ 
gle) the X-ray emission could possibly come from 
the outer magnetosphere, with an antenna pattern 
(and a beaming factor) possibly similar to the 7 - 
ray one. This could explain the narrow F.^/Fjc dis¬ 
tribution for RQ pulsars in spite of their expected 
large dispersion in 7 -ray beaming factors. Focus¬ 
ing on the RL pulsars, about half of the sample has 
Log(F.y/Fx) ^ 1 . 8 , much lower than the average 
Log(F.y/Fx) ^ 3.4 for the RQ pulsars. This could 
suggest the existence of a luminous X-ray emission 
component that can only be seen for small mag¬ 
netic impact angles, and thus possibly centred on 


the magnetic axis. Marelli et al. (2014b) explained 


the peculiar multi-wavelength behavior of the RQ 
pulsar J1813—1246 by modeling its X-ray emission 
using a polar cap model ( Dyks et al.|2004 ): a low- 
altitude cone beam with peak emission just inside 
the polar cap rim. The radio beam, centred on 
the X-ray cone, is missed by just a few degrees: 
a slightly different line of sight would have made 
J1813—1246 a bright radio pulsar. Such an X- 
ray polar cap emission component should display 
a large variability in luminosity and/or beaming 
factor within the RL pulsar sample in order to ac¬ 
count for the large dispersion in F.y/Fx. Such a 
variability could be driven by other pulsar proper¬ 
ties such as the spin-down-luminosity and/or the 
magnetic field configuration/intensity, as well as 
their evolution as a function of time. 

A major step forward in our understanding of 
the overall multi-wavelength emission geometry 
of pulsars could be obtained by performing si¬ 
multaneous fitting of their multi-wavelength light 
curves. Recently, Pierbattista et al. (2014) have 


jointly fit 7 -ray and radio light curves with simu- 











































lated 7 -ray and radio emission patterns. This al¬ 
lowed them to investigate some relations between 
observable characteristics and intrinsic pulsar pa¬ 
rameters, encouraging the creation and testing of 


new models (see e.g. Kalapotharakos 2012 Li et 
ahjlM^ . Also owing to the increasing number of 


RQ pulsars, the above analysis should be extended 
to include X-ray light curves, taking into account 
both thermal and non-thermal pulsed emission. 
Simultaneous fitting and phase-resolved spectral 
analysis of multi-wavelength emission of pulsars 
will allow us to test different emission models and 
to build new models for the overall emission, also 
better explaining the results presented in this pa¬ 
per. Unfortunately, of the 28 RQ LAT pulsars now 
detected in the X-rays, only seven have been ob¬ 
served for a sufficiently long integration time to de¬ 
tect X-ray pulsations. Similarly, of the 30 RL LAT 
pulsars with an X-ray counterpart, only 15 are 
known to pulsate in the X-ray band. Therefore, 
an important contribution to the theoretical anal¬ 
ysis would come from the detection of X-ray pul¬ 
sations for most of the brightest pulsars through 
deep follow-up XMM-Newton and Chandra obser¬ 
vations. While a number of much brighter RQ 
and RL pulsars are better suited for such a future 
study, deep observations of XMM-Newton , the 
best in-flight telescope to perform such a search, 
could allow us to detect pulsations from our eight 
pulsars. For instance, if we assume a Lorentzian 
light curve with a duty cycle of 0.1, with an entire 
orbit of XMM observation (130 ks) we would de¬ 
tect pulsations at ha in case of a pulsed fraction 
of ^40% for the brightest of our eight pulsars to 
100% for the weakest one. 
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Fig. 1. — Chandra and XMM-Newton images of the pulsars listed in Table 1, taken in the 0.3-10 keV energy range. For a 
better visualisation, we applied a Gaussian filter with a kernel radius of 3^^ In each panel, the best-fit X-ray position (90% 
confidence errors) is shown by the white circle and the LAT 7 -ray timing position (Icr errors) by the green ellipses. The latter 
are taken from the most recent compilation of the LAT 7 -ray Pulsar Timing Models. The position of PSR J1838—0537 is taken 
from lPTetsch et al.| l |2012a[ |, where the positional error is not reported, and is marked by the green cross. 
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Fig. 2.— Histogram of the logarithm of the 7 -to-X-ray flux ratio (F^/Fx) of LAT pulsars with high-confidence X-ray 
detections, as defined in Figure 3 of |Marelli et al?] | |2011[ |. Histograms for RL and RQ pulsars are shown in red and blue, 
respectively. The increment to the radio-quiet pulsar histogram for the eight pulsars that we detected in the X rays for the first 
time is shown in dashed blue and is added to the histogram of previously-known radio-quiet pulsars. For these pulsars, we used 
the results of the PL fits to compute the Fx values. The continuous and dashed lines report the best gaussian fit for the RL 
and RQ distributions, respectively. 


Name 

RA(J2000) 

h m s 

Dec(J2000) 

O ' " 

P 

ms 

P 

10-14 s s-i 

E 

10^4 ej.g g-i 

Tc 

ky 

D.,, 

kpc 

J0734-1559 

07 34 45.7 

-15 59 18.8 

155 

1.25 

13.2 

200 

1.3 

J1429-5911 

14 29 58.6 

-59 11 36.0 

116 

3.05 

77.4 

60 

1.7 

J1838-0537 

18 38 56.0 

-05 37 09.0 

146 

46.5 

593 

5 

1.8 

J1846-f0919 

18 46 25.9 

-H09 19 45.7 

226 

0.99 

3.4 

360 

1.4 

J1957-f5033 

19 57 38.4 

-h50 33 21.4 

375 

0.68 

0.5 

870 

0.8 

J2028-f3332 

20 28 19.9 

-^33 32 4.2 

177 

0.49 

3.5 

570 

0.9 

J2030-f4415 

20 30 51.4 

+U 15 38.8 

227 

0.65 

2.2 

550 

0.8 

J2139-f4716 

21 39 55.9 

+47 16 13.0 

283 

0.18 

0.3 

2500 

0.8 


Table 1: Characteristics of the 7 -ray pulsars discussed in this work. Here, we list the name, 7 -ray timing 
coordinates, period (P), period derivative (P), energetics (P), characteristic age ((r), and pseudo-distance 
(D^), respectively. These values are taken froni |Abdo et al. (2013). Pulsars positions are taken from the last 
LAT 7 -ray Pulsar Timing Models, where possible, and from Pletsch et al. (2012a I for PSR J1838—0537. 
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Fig. 3. — Histogram of the best fitted 7 -to-X-ray photon index ratio (P-y/Fx) of LAT pulsars with high-confidence X-ray 
detections, as defined in Figure 3 of |Marelli et al^ | |2011[ |. Histograms for RL and RQ pulsars are shown in red and blue, 
respectively. The pulsars for which the photon index has not been fitted are not included. The continuous and dashed lines 
report the best gaussian fit for the RL and RQ distributions, respectively. 
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